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Abstract 



PhI We develop a relativistic quark model for pion structure, which incorporates the 

p ' non-trivial structure of the vacuum of Quantum Chromodynamics as modelled by 

^ ' instantons. Pions are boundstates of quarks and the strong quark-pion vertex is 

determined from an instanton induced effective lagrangian. The interaction of the 
constituents of the pion with the external electromagnetic field is introduced in 
^ I gauge invariant form. The parameters of the model, i.e., effective instanton radius 

c3 I and constituent quark masses, are obtained from the vacuum expectation values of 

the lowest dimensional quark and gluon operators and the low-energy observables 
of the pion. We apply the formalism to the calculation of the pion form factor by 
means of the isovector nonforward parton distributions and find agreement with the 
experimental data. 



Pacs: ll.15.Kc, 12.38.Lg,12.39.Ki,14.65.Bt 
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1 Introduction 



The quark and gluon distribution functions play an important role in the exploration 
of hadron structure. Perturbative quantum chromo dynamics (pQCD) allows one to 
obtain the evolution of the distribution functions with using operator product ex- 
pansion and renormalization group methods from some starting input, which cannot 
be calculated from the first principles. We have been therefore motivated to develop 
an effective approach which, based on fundamental principles of QCD, serves to 
estimate the observables associated with the dynamics of large distances. 

Recently much attention has been paid to the off-diagonal quark (and gluon) 
distribution functions parametrizing asymmetrical hadron matrix elements of the 
light-cone quark-gluon operators. The off-diagonal distribution functions are a gen- 
eralization of the conventional distribution functions, which serve as a link between 
the hadronic structure functions and their form factors. 

In here, we calculate the isovector off-diagonal leading-twist quark distribution 
functions for the pion within a novel approach, which is an extension of the well- 
known instanton model incorporating bosonization. An effective quark-hadron 
lagrangian of nonlocal type, where the nonlocality is generated by instantons, un- 
derlies our approach. We implement gauge invariance with respect to external elec- 
tromagnetic field in the effective lagrangian and the Green functions by introducing 
explicitly the Schwinger phase factor in the definition of quark fields. The parame- 
ters of the instanton vacuum, the instanton radius and the constituent quark mass, 
are related to vacuum expectation values of the lowest dimension quark-gluon opera- 
tors. The coupling constants are calculated from the compositeness condition which 
is formulated in our case as the vanishing of the renormalization constant of the 
meson-field. The quark distribution functions obtained correspond to a hadronic 
scale (point of normalization) given by the typical relative quark momenta that 
flows through the instanton vertex and is of order /io ~ P ~ p^^ ~ 0.5 -i- 1 GeV . We 
ought to stress that at some point in our development we restrict the calculation to 
the isovector part of distribution which is the only one required by the scrutinized 
observable [} 

We proceed in the next section to describe the nonlocal four-quark model induced 
by instanton exchange in the effective single instanton approximation. 

Section 3 is devoted to fix the parameters of the model from the experimental 
values of the vacuum condensates and the properties of the pion. In section 4 we 
calculate the corresponding nonforward parton distributions associated with deeply 
virtual Compton scattering and determine the pion form factor. In an appendix we 
explicitly show the contributions from the different diagrams. 



A more complete analysis in a different model can be found in ref. 
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2 Model 



The investigation of the ground state properties of QCD has important consequences 
for the description of the hadronic processes within the theory. At present, there 
exists at least three main approaches for studying the structure of the vacuum: i) the 
quasi-classical approximation; ii) a formalism based on Wilson's operator product 
expansion ; iii) lattice QCD. As it is well-known, the main problem of the quasi- 
classical method is that, for the 3 + 1- dimensional Yang-Mills theory, the instanton 
gas approximation is not able to provide the correct vacuum- vacuum amplitudes in 
the background of large QCD-vacuum fluctuations. In a scale-invariant theory, like 
QCD, the vacuum fluctuations can have large sizes leading to an infinite instanton- 
anti-instanton medium density, this is the so called infrared divergence problem, 
which hinders the solution. 

At the present, one way for stabilizing the pseudoparticle medium is the instanton 
liquid model ( see, for example, []T], In the framework of instanton liquid model, 
the medium of / — J is quite rarefied, i. e., the ratio of the average size to the average 
distance between pseudoparticles being approximately equal to 1/3. Besides, in the 
limit {Nc oo) the size of all (anti-)instantons is taken to be equal to some average 
size p. The instanton model with fixed average size of instantons and instanton 
density is the effective single instanton approximation of the QCD vacuum . 

Within the single instanton approximation the chirally invariant nonlocal four- 
fermion action has the following form 

5'4g = ^Gl J d^Xo J Y[d'^XnKi{xi, X2, X3, Xi) 

1^ [QniXo - xuXo)T,Ql{Xo, Xq + x^)] [Qn{Xo - X2, Xo)T,Ql{X, X + x^)_ 

+iR^L)} (1) 

where Xq is the position of the center of the instanton, and the matrix combinations 

Fj X Ti are 

(r^'^r^), ^(^^^_^^ (rV,,^rV,.), (2) 

t"" = (lyir) are the matrices characterizing flavor space and Xc = 3 is the number 
of colors. To ensure the gauge invariance of the nonlocal action with respect to 
an external electromagnetic Afj,{z) field, we define, following 0, H, the chiral quark 
fields with the Schwinger phase factor 

QR{L){x,y) = ^ ^^^ Q{x,y) 

Q{x,y) = exp -ieQ / dz^A^{z) q{y) 



X 



Q = diag(|,-^,-^~). (3) 
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In the local limit the action (|l|) describes quarks interacting through the 't Hooft 
vertex |§| and it is invariant under global axial {q{x) — > exp (i75r ■ 6')g(x)) and 
vector {q{x) — >■ exp {ir ■ 6)q{x)) transformations, and it anomalously violates the 
Ua{^) symmetry {q{x) exp {i'-f56)q{x)). In what follows we neglect the current 
quark mass and restrict ourselves only to the nonstrange quark sector. In the 
kernel of the four-quark interaction Kj{xi, X2, x^, X4) has the form 

[Xi, X2, X^, Xi) = f{xi)fix2)f{xs)fix,), (4) 

where f{xi) are related to the profile function of the quark zero modes (see be- 
low). Morever, this kernel characterizes the range of the nonlocality induced by 
quark-(anti-)quark instanton interaction, which is determined by the average size 
of instantons. The latter represents a natural cutoff parameter of the effective low 
energy theory. 

The 4-fermion nonlocal action is linearized by introducing some auxiliary field 
interpreted as pion field. Since this bosonization procedure has been elucidated in 
the literature (see, for details, ^^), let us immediately turn to the linearized effective 
quark-pion interaction lagrangian in gauge-invariant form which is given by 0, 

C%{x)=g^ggJ «2/(ei)/(6Ma: + ei,x)275r-7r(x)Q(x,a;-6). (5) 

The physical pion-quark coupling constant (^^gg is determined from a so-called cou- 
phng constant condition []Tl] 

Z^ = l+ gl^^U'iml) = 0, (6) 

where Z^^ is the renormalization constant of pion field. The physical consequence 
of the condition = is that the pion field is always found in a dressed state. 
Besides, the coupling constant condition is a strong restriction because it implements 
the dynamics responsible for the formation of the bound state. 

Summarizing, lagrangian (^) and the coupling constant condition @ characterize 
our approach. 



3 Model parameters 

We briefly discuss the model parameters (for details, see ref. 10). Due to the effect 
of spontaneous breaking of the chiral symmetry, the momentum dependent quark 
mass Mq{k) is dynamically generated. It obeys the well-known gap equation ^ 

^In a similar manner we can write the interaction lagrangian for other mesons. 
Here and below, all Feynman diagrams are calculated in Euclidean space (fc^ — —k'^) where the 
instanton induced form factor is well defined, and where it rapidly decreases so that no ultraviolet 
divergences arise. At the very end we simply rotate back to Minkowski space. One can verify that 
the numerical dependence of the results on the pion mass and the current quark mass is negligible 
and can be safely ignored in the following considerations. 
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given by 

d% P{k)M,{k) _ n, 



/ 



(27r)4fc2 + M2(fc) 4N^Mg' 

where Mq{k), the momentum dependent quark mass, is defined by the dressed quark 
propagator: 

Sp\p)=^-M,Q{p), (8) 

and, ric = M^/Gj is the density of the instanton medium {Gj is the four-fermion 
interaction constant in (|l|) and Mg = Mg{0)). The solution of eq.(|^) has the form 

Mg{k) = MgQ{k), (9) 

where Q{k) = P{k) is the nonperturbative part of the quark propagator which in 
the axial gauge is given by 

Q{P) = ^^^Jd^xexpi-tp-x)Qix'), g(0) = l; (p = \p\) , (10) 

and 

Q{x') = (:g(0)g(a;):)/(:g(0)g(0):), Q(0) = 1 . (11) 



In the axial gauge Eq. ([ll|) coincides with the gauge-invariant correlator in which 
the quark fields are connected by the Schwinger phase factor. In practical calculation 
we shall use the approximation for (|l^, |TlD as in ref . |jl2 



f{k) ^ 27rp,[^2.25exp{-p,|fc|} - 1.25 exp {-3p,|A;|} j . (12) 

Given the value of the dynamical mass, one can obtain the value of the quark 
condensate, 

^ ^t.SM. - y) ^ -,N^j^^^^^M_, (13) 



and the average quark virtuality in the vacuum |]14|-|115 



2 ^ (: qPS :) ^ _AN^ f d^k M,{k) 

'~ (:gg:) {qq) J {2tiY k^ + M^ik)' ^ ' 

The ratio r] = Mg{Xg) / characterizes the diluteness of the instanton liquid vacuum. 

Here, we have to emphasize that in the axial gauge the gauge invariant quantity 
Mg{k) is expressed in terms of the gauge invariant nonlocal codensate ([TI|). Thus, 
the quark condensate and the average quark virtuality in the vacuum are defined in 
the gauge invariant manner. 
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For the moment, it is instructive to consider equations (0) - (|1^) neglecting 
M^{k), compared to fc^, in the denominator of the integrands. This approximation 
is justified in the dilute liquid regime where (fc^) ~ » M^A^). Using the 
explicit form in the momentum representation of the zero mode in the axial gauge 
g from dlD, we have H 



A 



Pc 



By inverting these relations, we are able to express the parameters of the instan- 
ton vacuum model in terms of the fundamental parameters of the QCD vacuum 



Pc 



TT 



(15) 



/^From these relations the model parameters, pc and Mg, can be estimated in terms 
of the quark condensate, (gg) ~ —(230 MeV)^ and the average quark virtuality. 
The latter was estimated using QCD sum rules, A^ = 0.4 ±0.2 GeV^ 0, and lattice 
QCD (LQCD) calculations, which yield A^ = 0.55 ± 0.05 GeV^ Q. These values 
lead to Pc ~ 2GeV^^ and Mq ^ O.SGeV. The combined analysis of the vacuum 
and pion properties confirms such estimates [112 1 and gives for the allowable range 
of parameters, 1.5GeV~^ < pc 
that the diluteness condition rj ■ 
window. 



< 2GeV-^ 0.22GeV < < 0.26GeV. We note 
= M^(A^)/Ag -C 1 is well satisfied within the whole 



4 Deeply virtual Compton scattering and nonfor- 
ward parton distributions 

Let us consider the deeply virtual Compton scattering (DVCS) within our model. 
The interest of such investigation is that it defines a procedure to calculate the 
nonforward parton distributions, which correspond to a generalization of the both, 
the parton distributions and the hadronic form factors. 

During virtual Compton scattering, a pion with momentum p absorbs the virtual 
photon of momentum q, producing an outgoing real photon of momentum q' = q + r 
and a recoiling pion with momentum p' = p — r. We shall concentrate on the deeply 
virtual kinematic region of g, i.e., the Bjorken limit: = —q^ — > oo, p ■ q —>■ oo 
and Q"^/ {p ■ q), t = r"^ finite. 

The expression for the Compton scattering amplitude of the photon off the pion, 
with momenta p, q in the initial state and p', q' in the final state, is given by 

T^,y{p, q) = j diidi2dizdii exp {iqii - iq'^2 + ip^s - ip'U) 

6^S 



)o. (16) 
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The S-matrix has the standard form 



S = Texp (^i J dx [Ceff{x) + Cem{,x)]^ . (17) 

This amphtude corresponds to a set of diagrams drawn in figs. 1. Fig. 1.1 shows 
the contribution of the a meson to the process we are calculating, which is small, 
as can be shown numerically and by counting^. Moreover it can also be shown 
that the diagrams labelled (c) and {d) are suppressed in the Bjorken limit, therefore 
we shall only pay attention to the diagrams (a) and (6). Moreover, diagram (a) 
gives the biggest contribution. The general method for obtaining of the off-diagonal 
distribution function can be found in |TP|, W^. Following these methods, we describe 



in detail only the contribution of diagram (a) to the scattering amplitude in the 
momentum representation 

tr f75S(k - p)75S(k)7^S(k + q)7,S(k - r)) . (18) 



In order to obtain the distribution functions, it is convenient to define a special 
coordinate system (see, for instance |18[) in which light-like vectors p and q' are 



collinear. Then other vectors are expanded in terms of p, n = q' / {p ■ q') and the 
transverse vectors 

= (1 - OP/. - + P'ft^ = (pf, + Ui^ + r^, 
% = -CPf^ + (P ■ ?)^M + C = r\\ ■ n = 1 - p' ■ n, 

p-q = p-q' + ^, (19) 

where ( is the asymmetry parameter which in the Bjorken limit coincides with the 
Bjorken variable xb for deep-inelastic scattering 

xs = -^ = ^ 
2pq 2pq 

Now we insert the given decompositon of 4-vectors in eq. (|TBp and carry out 
the integration over k, keeping only terms which are not suppressed in the Bjorken 
limit. Then, using 

1 

J dX6{X - kn) = 1, (20) 



'*The direct a exchange diagram, which would be leading order in 1/N, does not contribute to 
the isovector OFPD. 

^ To define the parton distribution we formally turn to the Minkowski signature. 
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we obtain for the contribution of diagram (a) th the Compton scattering amphtude 
T^u{Q. q', r) = l {P,nu + n,Pu - 9,.) j dX^^T^{X, t) (21) 



where 



ti{l^S{k - p)^5S{k)h ■ ^S{k - r)^5(X - kn) + (X ^ X). (22) 

In a similar manner we obtain the contribution of diagram (6) which is given by 

1 r 1 

where 



.2 , j4; 1 



-^C(^,t) = ^'f^J^J{k)f{k + p')Jdrri{k-rrr + r^ril-r)) 



tr(^75^(A:)75^(A; + p')) - 

(^6{X ~kn) + 6{X -C + kn)j +{X ^ X). (24) 

The variable X is the total fraction of the initial hadron momentum p carried 
by the active quark, and satisfies the kinematical constraint < X < 1. 

In order to calculate the structure integrals in eqs. ( P^ and ( pij ) we use the 
a— representation for all propagators, the Laplace transform for the nonlocal vertex 
function f{k) and the integral representation for 5— function. Note that, in the 
/c— integrals which appear in the above equations, we use the Euclidean signature, 
but at the final stage we move back to Minkowski space for the external momenta. 
In the case of t = we consider a family of asymmetric distribution functions JF^(X) 
whose shapes change with (. The ( -dependence of the asymmetric function is the 
main difference with respect to the double distribution function F{x, y) recently 
introduced in |jl9|. The asymmetric function JF^(X) satisfies the sum-rules, 

1 

J dXJ^^{X) = 1. (25) 



This distribution function has the following properties [|l^]: When the total mo- 



mentum fraction X of the initial hadron momentum p is larger than the fraction 
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( of the momentum transfer r, the function JF^(X) can be treated as a generahza- 
tion of the usual distribution function; when X < (, the asymmetric function can 
be treated hke a distribution amphtude $ of the gg-states, in our case the pion, 
with momentum (p; besides, for C = 1, the asymmetric quark distribution function 
^^(X) becomes the pion wave function (see, for example, pO|). 

The asymmetric quark distribution functions J-'^ and JF^? are shown in fig. 2. 

Finally, the integration of Eqs. (^) over the total momentum fraction X 

yields the following sum rule: 

1 

JdXJ^^{X,t)=F^{t), (26) 



where -F^(t) is the hadronic form factor normalized to unity. The behavior of the 
form factor -F^(t) is shown in fig. 3. Our numerical results are in agreement with 
experimental data and also coincide with results given by other approaches (see, for 
example, |]2l[| ). 



5 Concluding remarks 

We have used the instanton liquid model of the QCD vacuum to obtain an effective 
quark-hadron lagrangian by bosonizing the 4-fermion interaction generated by the 
instantons. The procedure imposes a coupling constant condition which allows us 
to calculate the physical value of the quark-hadron interaction constant. 

The parameters of the model, related to the properties of the vacuum structure, 
have been fixed by using the < qq > condensate and the average quark virtuality of 
the vacuum. Furthermore, the pion properties confirm our estimates. 

We have considered within the model deeply virtual Compton scattering a pro- 
cedure which allows the calculation of nonforward parton distributions. We have 
calculated the asymmetric distribution functions at different asymmetrical parame- 
ters {(). Finally the integration over the total momentum fraction on the nonforward 
parton distribution yields the pion form factor. Our results agree with those of other 



authors and with experiments ||22 



This work represents step forward in the direction of deriving a model for hadron 
structure directly from QCD with well established assumptions. The procedure 
developed can be generalized to other mesons and to baryons, and the work in this 
direction is in progress. 
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Appendix A 



1. The contribution arising from the basic diagrams (see, fig. la): 



where 



1 °° 





-exp -t{l-X)- 



tti + Q;34(l-C) V Qii + a34(l-C)// 

oo 

dai..da5fi{ai)f2{a2)fi{a3) ^^ ^ ^)to35^6(«3) 







(q;i4 + q;35(1-C))' 



exp ( -M2(a45 + «3) - t(l - X)- 



Q:i4 + a;35(l - C), 

2. The contribution arising from the tadpole diagrams (see, fig. lb): 

2 - C 47r2 
where 



X^''^\Q;,i) = J dai..da3a3fi{ai)f2{a2)fi{a3) J dr ^ _ 



asCr - a2{l - C) 



e(a.)exp f a.M^ - tra, + '^^ + "l^^ ~ V 

\ Q;123 + ^2 / / 

T}'^\a,t) = ! dai..da2.aJ\{ai)j2{Q.2)J\{a3) I dr .^'^'^ . 







0(Q;3)exp dsM — tra^ 
V V '^123 + Q;3 

9(d4)exp q;4M^ — tra^ 



ai23 + 014 
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The notation is: 



Q^ni..ni Oiji^ -|- .. -|- Cy^ni 

1-x x-c 

Oil = "^2^ ^^^^C 

X X-C 

X X-C 
^ "'^^ 1 _ X "'^^ 1 - X ~ 

"1 = "3 ^ «2 «!, 

C(i - ^) - ^ 1 - ^ 

Q;2 = -0(3 C- X 0^2 ^ _ ~ '^l' 

C(i-t)-x c-^ 

as = -0(3 1_ X ^^ 1- X ~ 

X-Ct X 
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Figure 1: The diagrams contributing to the Compton Scattering amphtude: (a) is 
the basic type and (6) — (d) are tadpoles. 




Figure 2: The sigma exchange diagrams contributing to the Compton Scattering 
amphtude. 
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Figure 3: The asymmetric distribution functions at different values of an asymmetry 
parameter Tlie solid line shows the contribution to the distribution function of 
both the basic and the tadpole diagrams; the dashed-line represents the contribution 
from the tadpoles alone. 
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Figure 4: The pion elctromagnetic form factor. The sohd hne shows the contribution 
from both the basic and tadpole diagrams, the dashed-hne again represents the 
contribution from the tadpoles only. 
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